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Global  or  relaxation  formulations  for  the  re¬ 
duced  form  of  the  Navier-Stokes  equations,  frequent¬ 
ly  referred  to  as  parabolized  Navier-Stokes  (PNS) , 
are  presented.  Difference  procedures  and  relaxation 
solutions  for  the  (u.v.p)  system  are  presented. 

The  continuity  equation  is  satisfied  exactly  at 
each  grid  point  and  a  poisson  pressure  equation  is 
not  required  explicitly.  The  development  of  a  se¬ 
cond  composite  (U,$,G)  velocity  relaxation  procedure 
for  the  primitive  variable  equations  is  also  dis¬ 
cussed.  For  the  (u.v.p)  system,  several  model 
problems,  e.g. ,  finite  flat  plate,  trough ,  boattail 
and  airfoil,  are  considered.  Strong  pressure  inter¬ 
action  is  evident  in  each  case  and  axial  flow  se¬ 
paration  occurs  for  several  of  the  problems.  The 
questions  of  accuracy,  stability, convergence  rata, 
and  implied  difference  forms  of  the  pressure  and 
vorticity  equations  are  addressed. 

1.  Introduction 

Conventional  methods  for  the  numerical  solu¬ 
tion  of  the  primitive  variable  form  of  the  incom¬ 
pressible  (elliptic)  Navier-Stokes  or  ("semi- 
elliptic")  parabolized  Navier-Stokes  (PNS)  equa¬ 
tions  are  such  that  the  velocity  components,  u,v, 
are  determined  from  the  momentum  equations,  and  the 
pressure  p  is  obtained  from  the  differential 
poisson  equation  derived  from  the  momentum  aqua¬ 
tions.  The  equation  of  continuity  is  not  evaluated 
explicitly  but  is  satisfied  indirectly  through  the 
poisson  equation  and  pressure  boundary  conditions. 
Since  this  procedure  differs  markedly  from  most 
inviscid,  boundary  layer  and  triple  deck  formula¬ 
tions,  an  alternate  development  that  more  closely 
follows  these  asymptotic  theories  is  considered 
here  for  the  evaluation  of  viscous  interacting 
flows  at  large  Reynolds  numbers. 


The  objective  of  the  present  development  is 
the  solution  of  the  PNS  system  by  direct  application 
of  the  momentum  and  (first-order)  continuity  equa¬ 
tions.  The  formulation  does  not  require  the 
second-order  differential  form  of  the  poisson 
pressure  solver.  A  global  line  relaxation  procedure 
is  developed  for  (u.v.p)  or  a  composite  [6,  7] 

(U, <P ,G)  system.  For  non-separated  flows  only  p 
or  4  are  stored  during  the  relaxation  process.  For 
separated  flows,  (u,v)  or  (U)  is  required  only  in 
regions  of  reversed  flows.  This  significantly 
reduces  computer  storage  requirements. 


Governing 


We  consider  here  the  reduced  set  of  PNS  aqua¬ 
tions  (1,  2)  written  in  two-dimensional  or  axi- 
symmetric  body  fitted  conformal  coordinates.  The 
equations  in  general  orthogonal  coordinates  are 
given  in  [3] .  As  discussed  for  cartesian  coordinates 
in  [2,  3],  for  incospressible  flow  a  consistent  PNS 
approximation  allows  for  the  neglect  of  all  axial 
(O  diffusion  terms  as  well  as  all  diffusion  effects 
in  the  normal  (n)  momentum  equation.  Normal 
diffusion  can  be  included  in  the  n -women turn  equation) 
however,  for  consistency  these  terms  have  generally 
been  neglected.  Numerical  tests  with  and  without 
these  terms  have  confirmed  the  validity  of  this 
approximation  for  several  of  the  problems  consi¬ 
dered  herein. 


(i)  (u.v.p) ! 


continuity 


(hhjU) ^  +  (hh3v) n  ■  0 


C-momentum 


In  the  present  paper,  the  authors  continue  the 
line  of  thought  first  presented  for  the  PNS  system 
in  references  [1-3] .  The  analysis  is  developed  in 
greater  detail  here  and  in  [4,  SI.  The  questions 
of  global  stability  of  the  relaxation  procedure, 
the  resulting  difference  forma  of  the  pressure  and 
vorticity  equations,  accuracy  and  rate  of  con¬ 
vergence,  are  examined  more  critically.  Compari¬ 
sons  are  given  with  triple  deck  and  interacting 
boundary  layer  solutions  for  trailing  edge  and 
trough  configurations;  solutions  are  also  obtained 
for  boattail  and  airfoil  geometries.  The  effects 
of  strong  pressure  interaction  and/or  flow 
separation  are  evident  in  each  of  these  problems. 
The  majority  of  the  solutions  are  for  laminar  flow 
conditions;  however,  several  results  have  been 
obtained  with  the  Cebeci-Smith  two  layer  eddy 
viscosity  closure  model. 


(hh,u2) , *  (h2uv)  ♦  uvh.h  -  v2h,h. 
3  5  n  3  n  3  5 


n -momentum 


hhjp^  ♦  V-T/Re 


[h.(hu)  /h  ] 

3  n  n 


-  (hh,uv)  -  (hh,v2)  -uvh,+  u2h  h 

3?  3  n  5  3  n , 


The  cartesian  coordinates  5  *  5(x,y) , 
n  ■  n(x,y)  are  related  to  the  (x,y)  physical  coor¬ 
dinates  through  the  transformation  j«f(z)  or 
z«F(a)  where  a  -  5+in  and  z  -  x+iy.  The  metric  h 
is  defined  by 


This  research  was  supported  by  -he  Air  Force  Office  of  Scientific  Research  under  Grant  No,,  90-0047. 


88  04  ppr  r: 

diotrU-.-.t  :.o.:  i.; 


i  Ai-C 


1 


1/2  1/2 

h  -  |f(z)|  -  lx.*y.)  -  (x  +y„) 

t  %  n  n 


the  full  Navier-Stokes  equations  in  (6,  7).  Analysis 
of  the  PNS  equations  with  the  composite  system  is 
in  progress  and  results  with  this  formulation  shall 
be  presented  in  a  future  paper. 


The  metric  h.  «  y£,  where  £  »  0  for  two-dimensions 
and  £  •  1  for  axisyuostry .  The  metric  h  and  all 
derivatives  are  evaluated  with  second-order  differ¬ 
ence  formulas.  In  the  axial  and  normal  momentum 
equations  (lb.  lc) ,  the  metric  h  and  derivatives 
of  h  are  assumed  to  be  at  most  of  order  one.  For 
geometries  with  larger  curvature  the  complete 
expression  for  V-T  may  be  required  in  (lb)  and  the 
viscous  effects  in  (Id)  may  become  important.  The 
full  Navier-Stokes  equations  or  a  more  appropriate 
non-conformal  coordinate  mapping  may  then  be 
required,  see  (6,  7]. 

(ii)  (U,4.G) 

In  the  composite  velocity  development  described 
in  [6,  7],  an  "inviscid"  pseudo-potential  4 , 
"viscous"  velocity  a  and  "inviscid"  Bernoulli 
pressure  G,  replace  (u.v.p) ,  i.e., 

u  -  0(1  ♦  ),l/h  •  J  u( 

v  ■  *  /h 
n 

G  »  p/o  ♦  (u2  +  v2)/2  -  G 

e  • 


The  equations  became 
continuity 

(h  U4.) .  +  (h,4  )  +  (h,U) -  -  0  (2a) 

Jit  jnn  J  t 

^-momentum 

u  +  {(hh,  u2(U2-U>],+  (hh.u  v(u-l)]  >/h2h, 
t  je  £  j  •  n  3 

+  u  vh  (0-1) /h2  +  u  u  .  (0-l)/h 
e  n  e  at 

-  -  G./h  ♦  V-T/R^  (2b) 


n-momentum 


G  -  -  (0-1)  ( (u2/2)  -  u20h  /h}  (2c) 

n  e  n  a  n 

This  multiplicative  composite  velocity  development 
is  patterned  after  matched  asymptotic  viscous- 
inviscid  flow  theory.  For  inviscid  irrotational 
flows,  0=1;  (2a)  then  reduces  to  the  potential 
equation  and  (2b, c)  lead  to  the  Bernoulli  equation, 

G  «  0.  For  boundary  layer  problems,  u#, 

(or  4r,  4..)  and  G  are  specified  with  boundary  layer 

s  S'* 

edge  conditions  and  (2a, b)  combine  to  determine 
(0, v) .  Interacting  boundary  layer  theory  combines 
elements  of  both  limits,  so  that  (2a, b)  form  a 
coupled  system  for  (0,4),  with  G  prescribed.  The 
equations  (1)  or  (2)  contain  all  the  terms  appearing 
in  each  layer  of  the  triple  deck  structure  (31. 

The  present  investigation  concentrates  on  the 
(u,v,p)  formulation  described  by  equations  (1) . 
Composite  solutions  have  already  been  discussed  for 


3.  Difference  Equations 

In  the  previous  analyses  [1-3] ,  it  was  shown 
that  if  the  system  (1)  was  forward  marched  in  the 
boundary  layer  sense,  i.e.,  backward  differences 
are  applied  for  all  £  derivatives  in  non-separated 
regions,  the  elliptic  pressure  interaction  would 
not  be  properly  represented  and  therefore  the 
exponentially  growing  Lighthill  departure  solutions 
would  appear  for  step  sizes  A£  <  (A£)  .  .  For 

cartesian  coordinates,  from  [2,  3],  we  find  that 

(d£)  . _  «  —  y„,  where  y„  is  the  location  of  the 
nun  it  m  m 

outer  boundary  y  -  yM.  Only  for  yH  <<  1  are 
accurate  solutions  possible  with  forward  marching. 

If  global  relaxation  or  multiple  sweep  marching 
is  used,  i.e.,  all  £  derivatives  of  velocities  are 
backward  differenced  in  non-separated  regions,  but 
some  form  of  forward  differencing  is  applied  for 
P£  in  (la) ,  the  elliptic  pressure  interaction  is 
recovered  and  the  departure  free  limit  d£ > 

is  removed.  Solutions  can  then  be  obtained  for 
A£  •*  0,  see  [2,  3];  the  numerical  procedure  is 
consistent  and  any  dasired  degree  of  accuracy  can 
be  specified.  Finally,  in  order  to  circumvent  the 
pressure  singularity  at  separation,  the  pjr  term 
must  also  allow  for  a  local  as  well  as  a  spatial 
interaction.  For  example,  central  differencing 
fails  in  this  regard  and,  as  discussed  in  [2,  3], 
is  unstable  globally.  Forward  differencing  of  p^ 
satisfies  all  constraints  and  moreover  is  consis¬ 
tent  with  the  eigenvalue  analysis  of  Vigneron 
et  al.  [9]  which  shows  that  for  incompressible  flow 
(M  *  0) ,  there  should  not  be  any  forward  marched 
component  of  the  p-  term;  i.e.,  u  «  0  in  his 
analysis .  Forward ’differencing  and  global  relaxa¬ 
tion  was  first  applied  successfully  in  [1-3]  for 
several  model  incompressible  flow  problems.  The 
extension  to  compressible  flows  is  discussed  for  a 
conical  geometry  in  [2]  and  for  flows  with  axial 
flow  separation  and  strong  pressure  interaction 
in  [10] .  More  detailed  discussion  and  results 
are  given  in  [3,  4]. 

The  difference  scheme  used  in  [1-3]  was  devel¬ 
oped  from  the  following  discrete  grid: 


j-1  .  u,v,p 


u: 

i-l 

i 

i+1 

v: 

i-1 

i 

i+1 

p: 

i-l 

i 

i*l 

Figure  1:  Difference  Grid  I. 
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The  continuity,  x  end  y  momentum  equations  are 

centered  at  ©,  0  and  0,  respectively  and  the 

difference  equations  given  for  a  uniform  mesh  in 
cartesian  coordinates  are  as  follows. 

continuity,  centered  at  (i.j-*j): 


-  0  (3a) 

x-aomentum,  centered  at  (i,j): 


■/-momentum,  centered  at  (i,j->j): 


The  equations  are  shown  here  in  non-conservative 
form  *f  this  simplifies  the  subsequent  discussion 
of  the  linear  system,  in  fact,  for  most  calcula¬ 
tions  conservative  equations  were  considered. 

All  quantities  are  evaluated  at  the  ith  iteration 

level  except  p“+*»  which  is  evaluated  at  the 
previous  iteration.  For  separated  flows  the  con¬ 
vective  terms  are  upwinded  and  terms 

are  also  required  at  level  n-1.  The  difference 
equations  (3)  are  first-order  accurate ,  i.e., 

0(4x,  iy2). 

An  alternate  and  more  accurate  derivation  of 
the  equations  and  interpretation  of  "forward" 
differencing  for  p^  is  given  below  for  cartesian 
coordinates.  This  system  was  considered  initially 
for  inviscid  flows  [4]  and  resembles  a  slightly 
different  development  proposed  by  Israeli  (8] . 
Consider  the  staggered  difference  grid  as 
shown: 


j 

j-1 


0  ©©  w 


V 


u:  i-1  I 

v:  i  i+1 

p«  i  i+1 

Figure  2:  Difference  Grid  II 

The  appropriate  difference  equations,  also  shown  in 
non-conservative  form,  are  now  centered  at  somewhat 

different  ©  ,  0,0  locations, 
continuity,  centered  at  ©  is  the  same  as  (3a): 


(4b) 


y-momentum,  centered  at  0  : 
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In  this  formulation  full  second-order  accuracy  is 
achiaved.  Tha  primary  modification  of  tha  systaa 
(3)  is  tha  averaging  of  tha  y-derivative  tanas 
in  tha  momentum  aquations.  Tha  unknown  pressure  p^ 
is  also  shifted  one  point  to  the  left  of  that 
given  by  the  formulation  (3).  This  interpretation 
is  more  accurate  and  is  consistent  with  the 
character  of  the  interactive  solutions  as  will  be 
seen  for  tha  trailing  edge  problem  to  be  discussed 
in  a  following  section. 

A  third  system  of  equations,  which  also  pro¬ 
vides  second-order  accuracy  has  been  proposed  by 
Israeli  IS].  The  staggered  u,v  grid  is  then  of 
the  form: 


j+l 

• 

U,P 

V 

U> 

V 

j 

u,P 

©® 

U,p 

5-1 

© 

© 

U) 

UfP 

V 

u: 

i-1 

i 

w  term  in  vw  ;  however,  with  tha  definition  (5b) 

y  2  2 

tha  error  in  this  expression  is  0(Ax  ,  Ay  ) , 
so  that  sacond-ordar  accuracy  is  retained  in  tha 
difference  approximation  of  the  vorticity  transport 
equation,  sea  [4] .  Similar  results  can  be 
obtained  with  tha  scheme  proposed  by  Israeli  in 
Figure  3. 

If  (5b)  and  the  continuity  aquation  (3a)  are 
combined  to  eliminate  either  u  or  v  terms,  for 
u  «  0  we  recover  a  nine-point  second-order  accurate 

2  2 

difference  formula  for  either  7  u  ■  0  or  7  v  •  0. 

In  fact,  inviscid  irrotational  flows  can  be  solved 
numerically  with  (5b)  and  (3a),  in  lieu  of  the 
potential  equation  72$»o.  This  is  a  result  of 
backward  differencing  of  ux  in  the  continuity 
aquation  and,  as  seen  in  (5b) ,  forward  differencing 
of  vx  in  the  definition  of  vorticity.  One  boundary 
condition  is  satisfied  for  u  (left  boundary)  and 
one  for  v  (right  boundary) .  This  is  a  direct 
result  of  the  differencing  procedure  applied  for 
the  u,v,p  primitive  variable  system;  i.e. ,  backward 
differences  for  velocities  (in  non-separated 
regions)  and  ’forward"  pressure  differences  as 
interpreted  in  Figure  1  (first-order  accuracy)  or 
Figure  2  or  3  { second-order  accuracy) .  Further 
details  of  this  analysis  are  given  in  [4] . 


(ii)  Difference  aquation  for  pressure: 


Figure  3:  Difference  Grid  III. 


In  a  similar  manner,  the  effective  poisson 
difference  equation  for  the  pressure  can  be  obtained 

by  (x-o«)ifj-  («-*»)i.1,jMy-mom)i  j+1-  (y-«c-).f 

-  0.  This  equation  is  of  the  form  (see  [4]  for 
derivation ) 


The  equations  are  only  slightly  modified  from  those 
of  (4) ;  however,  the  y  boundary  condition  for 
v  must  be  treated  remewhat  differently. 

(i)  Definition  of  vorticity  and  vorticity 

transport  equation: 

If  the  nonlinear  coefficients  u  aiJ  are 
assumed  constant  and  the  pressure  is  eliminated 
from  (4a, 4b),  the  following  difference  equation 
is  obtained,  in  the  inviscid  limit,  for  the  vorti¬ 
city  transport  [4] : 


“i . j+l^i .  1-l**‘‘i-l.i+l~<*ll-l .  1-1. 


°*i.j*l-  ffPi"j  *pi.j>  *2*,pi.j”  +Cpi.j-! 
♦pW.j*pi-l,j  ’  fi,j(u’V) 


where  a  -  (Ax/A y) 

In  order  to  recover  the  more  conventional  line 
relaxation  form  of  the  difference  equation, 

Israeli  (5]  has  shown  that  if  a  source  term  S.  . 

» ,  J 

is  introduced  in  the  x-momentum  equation  (3b  or  4b), 
we  obtain 


«  0(Ax2,  Ay2) 


api , 5*1  '  2 <1+*)pi, 3  +  °pi , 5-1 +  Pi*l . j  +  pi-l, 5 


fi.^'V) 


(5b) 

This  corresponds  to  defining  w  at  the  location  u 

shown  in  Figure  2  and  centering  (5a)  at  location  ©. 
Both  expressions  (5a)  and  (Sb)  are  second-order 
accurate.  The  difference  equation  obtained  by 
replacing  w.  ,  terms  in  (5a)  with  tha  formula  (5b) 

*  •  J 

is  exactly  that  resulting  from  the  elimination 
of  the  pressure  in  (4a)  and  (4b)  except  for  the 


This  is  equivalent  co  introducing  an  iterative 
"time''  derivative  into  the  relaxation  process. 
Israeli  [5]  has  also  introduced  an  overrelaxation 
parameter  c  as  a  mechanism  to  improve  convergence 
rate.  Our  experience  has  not  led  to  a  marked 
improvement  in  convergence  for  the  nonlinear 
system  (4)  with  over relaxation.  Therefore,  all 
solutions  depicted  herein  are  for  the  first-order 
system  (5a)  or  the  second -order  system  (5b)  with 
c  »  1 .  Multi-grid  procedures  have  been  applied 
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for  convergence  acceleration  when  fine  meshes  are 
required,  see  (4). 

It  should  be  reiterated  here  that  although  a 
poisson-like  relaxation  scheme  can  be  inferred, 
boundary  conditions  are  prescribed  for  p  only  at 
the  outer  (y)  boundary  and  downstream  (x)  boundary; 
also,  the  pointwise  continuity  equation  is  satisfied 
exactly.  As  noted  in  the  introduction,  this  is  in 
sharp  contrast  with  conventional  solution  procedures 
that  use  the  poisson  form  of  the  pressure  equation 
with  full  Neumann  boundary  conditions  and  only 
indirectly  satisfy  the  pointwise  continuity 
equation. 


4.  Consistency  and  Convergence 


In  order  to  complete  the  analysis  of  the 
relaxation  technique  discussed  herein,  two  questions 
are  posed.  Does  the  prescribed  differencing 
procedure  capture  the  elliptic  pressure  interaction 
in  each  sweep  of  the  iteration  cycle,  and  what  are 
the  convergence  properties  of  the  global  relaxation 
process.  The  first  question  has  been  addressed 
in  [1-3]  and  as  discussed  in  section  3,  the  depar¬ 
ture  effect  is  circumvented  and  the  elliptic  inter¬ 
action  is  captured  with  the  forward  pressure 
differencing  of  (3b)  or  (4b).  as  shown  in  [4],  the 


introduction  of  the  source  terms  S 


provides 


increased  numerical  damping  and  enhances  the 
stability  properties  in  each  marching  step. 


The  second  question  is  considered  in  detail  in 
reference  [4] .  The  primary  conclusions  are  pre¬ 
sented  here.  Although  the  forward  pressure  dif¬ 
ferencing  eliminates  the  inconsistency  found  with 
single  sweep  procedures,  so  that  AC  can  be  made 
arbitrarily  small,  the  departure  limit  Ax  >  2y^/s 

appears  indirectly  as  a  factor  affecting  convergence 
of  the  global  relaxation  procedure. 

It  can  be  shown  from  a  linear  global  stability 
analysis  that  central  differencing  of  px  is  unstable 
[2-4] ,  but  with  forward  differencing  of  px,  as  in 
(3b),  the  iteration  technique  is  unconditionally 
stable.  The  maximum  eigenvalue  is  given  (3,  4]  as 


i-vJ  f' 


A  -  1 


where  c,  is  a  constant  of  order  one,  y  is  the  outer 

*  n 

y  boundary,  xu  is  the  outer  x  boundary,  and  N  is 

M  X 

the  number  of  x  grid  points;  i.e.,  N  «  x„/Ax. 

y^.  X  n 

Th«r*for«  for  — —  <<l  or  Ax  <<  — ,  X  •+*  l  and 

Yjl  x 

the  convergence  rate  will  deteriorate .  For 
(~)  <<  1  and  N^  fixed,  convergence  is  very  slow 
as  A-l  -  .  Although  the  departure  limit 

'Ax  " 

r—  >  1  is  no  longer  a  stability  limitation,  the 
condition  — —  <<  1  is  a  convergence  limit.  For 


Ax  1  i 

co*r**r  grid*,  o-g.,  —  >  rrr  #  convergence  is 

yM 

quite  rapid.  For  finer  grids,  SOR  or  multigrid 
acceleration  hae  bean  considered  [4,  5] . 

The  outer  pressure  boundary  condition  must  be 
prescribed  at  a  location  beyond  the  extent  of  the 
interaction  zone  of  influence,  e.g.,  triple  deck. 

The  pressure  boundary  condition  can  be  fixed  during 
each  sweep  of  the  global  procedure.  This  value  will 
remain  unchanged  if  the  outer  boundary  is  suffi¬ 
ciently  far  from,  and  unaffected  by,  the  viscous 
interaction;  alternatively,  the  pressure  boundary 
condition  can  be  updated  prior  to  each  sweep  in 
order  to  account  for  viscous  displacement  effects. 

In  conformal  body  fitted  or  streaadine  coordinates, 
this  should  be  unnecessary.  Unlike  interactive 
boundary  layer  theory,  where  the  outer  pressure 
boundary  value  requires  a  local  interactive  treat¬ 
ment  in  order  to  circumvent  the  separation  point 
singularity,  a  fixed  outer  pressure  condition  is 
acceptable  with  the  PNS  formulation.  The  normal 
momentum  equation  reflects  the  outer  in viscid  or 
interactive  behavior  and  the  separation  singularity 
is  automatically  suppressed. 

To  simarize,  if  px  is  treated  explicitly  or 
central  differenced,  the  global  procedure  is 
unstable;  if  px  is  backward  differenced,  departure 
solutions  appear  for  Ax  ■*  0.  with 

"forward"  differencing  all  iterative  procedures  ue 
stable  and  the  global  marching  problem  is  well- 
posed. 


5.  Boundary  Conditions 

For  the  finite-difference  grids  of  Figures  1 
or  2  the  appropriate  boundary  conditions  are 
specified  as  follows  in  the  transformed  body  fitted 
coordinate  system: 

At  a  surface  y  -  0  ( j  »  0) ,  the  velocities 
u  ■  v  »  0;  at  a  synmetry  line  y  «  0  ( j  -  0) ,  the 
velocities  satisfy  u^  »  v  -  o. 

At  the  outer  boundary  y  •  yM  ( j  -  M) ,  where  yM 

lies  outside  of  the  extent  of  the  interaction 
zone,  e.g.,  triple  deck,  p  ■  p^,  u  «  u^  and  a 
boundary  condition  on  v  is  not  required. 

At  the  inflow  boundary  x  »  0  (i  ■  1) , 
u  »  u(0,y) ,  and  vx(0,y)  »  0.  For  invlscid 

regions,  where  u  ■  u^,  vx  “  0  is  a  zero  vorticity 
condition  and  for  viscous  regions  v^  »  0  is 

equivalent  to  a  boundary  layer  approximation. 

The  velocity  v  should  not  be  specified  at  the 
Inflow.  This  leads  to  inviscid  vorticity  produc¬ 
tion  and  has  a  destabilizing  effect  on  the  global 
iteration  procedure.  The  flow  pressure  is  not 
prescribed  and  with  the  formulation  of  Figure  2 
is  unknown  and  a  result  of  the  calculation 
procedure. 

Finally,  the  only  boundary  condition  required 
at  the  outflow  is  the  pressure  or  equivalent 
pressure  gradient.  This  of  course  reflects  the 
elliptic  pressure  interaction. 
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6.  Solution; 


Five  nodal  problams  hava  sarvad  u  Cast  casas 
foe  tha  global  PNS  formulation  dascribad  harain. 

For  aach  geometry  thara  ia  a  region  of  strong  pras- 
sura  interaction  and  in  savaral  caaas  axial  flow 
saparation  occurs.  Tha  tast  problams  include 
(1)  the  trailing  edge  of  a  flat  plate,  (2)  tha 
Carter-Wornom  [121  trough,  (3)  a  boattail  configu¬ 
ration,  (4)  a  NACA  0012  airfoil  at  zero  incidence 
(laminar) ,  and  (S)  the  NACA  0012  airfoil  at  zero 
incidence  (turbulent) . 

(i)  Trailing  Edge 

Solutions  for  the  trailing  edga  geometry  are 
given  in  figures  (4a, 4b).  The  agreement  with  tha 
interacting  boundary  layer  results  Of  [11]  are 
quite  good.  The  finest  grid  includes  (161x121) 
mesh  points  for  (x.y) ,  respectively.  The  coarsest 
grid  was  (41x121)  and  full  convergence  required 
only  several  global  iterations.  If  the  calculation 
was  run  on  the  finest  grid  alone,  convergence  was 
still  not  achieved  after  several  hundred  iterations. 
With  a  multigrid  technique  [4],  full  convergence 
to  0(10~4)  for  the  maximum  error  in  successive 
iterations  was  achieved  in  approximately  ten  to 
fifteen  global  iterations.  Tha  outer  boundary  yM 
was  chosen  to  lie  outside  the  triple  deck  extant. 

If  yM  violated  this  condition,  the  calculation 

diverged.  The  calculation  was  relatively  insensi¬ 
tive  to  yH  when  this  condition  was  satisfied. 

The  solutions  for  pressure  and  skin  friction, 
both  defined  with  triple  deck,  normalization  [111 

are  shown  for  Re  -  105.  It  is  significant  that 
the  skin  friction  (velocity  profile)  is  relatively 
insensitive  to  the  grid  and  appears  to  be  quite 
acceptable  even  on  seme  of  the  coarser  mashes. 

On  the  other  hand,  the  pressure  is  extremely  grid 
sensitive  and  requires  the  finest  mesh  in  order 
to  accurately  represent  the  triple  deck  interaction. 
With  the  difference  grid  of  Figure  1,  tha  minimum 
pressure  occurs  one  grid  point  downstream  of  the 
trailing  edge.  With  the  grid  of  Figure  2  this 
pressure  value  is  correctly  obtained  at  the  trail¬ 
ing  edge. 

( ii )  Trough 

The  solutions  for  tha  trough  geometry 
tyw(x)  -  e  sech  4<x- 2.5)1,  (0  <  x  <  -)  are  shown  in 

Figures  (5a,  5b).  values  of  t  *  -0.015  [3,  41  and 
€  »  -0.03  were  considered.  Only  the  latter  results 
are  presented  here.  Solutions  were  obtained  for 

Reynolds  msobers  up  to  Re  *  3.6  x  10*.  Again,  the 
agreement  with  the  interacting  boundary  layer 
solutions  is  quite  good.  The  insensitivity  of  Cf 
and  the  sensitivity  of  p  to  the  grid  is  also 
evident  for  this  example.  As  the  Reynolds  number 
was  increased,  more  smoothing  was  required  on  the 
coarser  grids  in  order  to  achieve  convergence  to 
the  prescribed  tolerance.  The  outer  undisturbed 
pressure  boundary  condition  was  held  fixed  at 
y  «  y  throughout  the  computation.  There  were  no 

difficulties  at  separation  or  reattachment  points. 

As  with  the  trailing  edge  problem,  full  convergence 
with  the  multi-grid  iteration  procedure  was 
achieved  in  ten  to  fifteen  global  iterations.  For 
the  finest  grid  (241x121)  mesh  points  were 
evaluated. 


1.10'  T.50 


Figure  4a.  Trailing  Edge  Pressure  Distributions 
PNS  Solver. 


Figure  4b.  Trailing  Edge  Skin  Friction  Solutions 
PNS  Solver. 
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(iii)  Boattail 


ni 


Figure  5a.  Trough  Skin-Friction  Solutions 
PNS  Solver:  E--0.03. 
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Laminar  flow  solutions  for  the  boattail  geo¬ 
metry  of  Figure  6  are  shown  in  Figures  (7a,  7b,  7c). 
The  grid  is  generated  with  the  Schwarz-Christoffel 
mapping  routine  of  Davis  [131 .  These  results  are 
in  good  agreement  with  full  Navier-Stokes  solutions 
obtained  with  the  composite  (u,$,G)  equations  as 
reported  in  [3,  7].  For  Re  ”  6000,  based  on 
maximum  radius,  with  a  juncture  angle  of  12  degrees 
a  sizable  separation  bubble  is  obtained.  All 
velocities  are  stored  in  the  recirculation  region. 
The  relaxation  process  is  slower  than  for  the 
trailing  edge  or  trough  geometries!  however,  con¬ 
vergence  to  10~4  for  the  maximum  error  in  pressure 
is  obtained  in  approximately  35  iterations.  As 
the  Reynolds  number  or  corner  angle  is  increased, 
the  rate  of  convergence  decreases  and  the  multi¬ 
grid  procedure  also  deteriorates.  Further  analysis 
of  this  behavior  is  required.  Some  improvement  has 
been  observed  with  the  source  correction  of  (5b) . 

Solutions  have  also  been  obtained  for  turbu¬ 
lent  flow  conditions.  The  Cebeci-Smith  two  layer 
viscosity  model  has  been  applied  to  close  the 
system.  Although  this  may  not  be  an  accurate 
approximation  in  the  recirculation  region,  it 
does  serve  to  give  a  qualitative  picture  of  the 
flow.  A  Reynolds  number  of  Re  *  5  x  10s  based  on 
maxi  mum  body  radius  has  been  specified.  The 
effective  turbulent  Reynolds  number  is  of  course 
much  lower  and  the  separation  region  is  considerably 
smaller  than  that  obtained  for  the  laminar  flow 
at  Re  «  6000.  These  results  are  discussed  in 
greater  detail  and  figures  are  presented  in  [4] . 

(iv)  Airfoils:  laminar  and  turbulent 

The  flow  over  NACA  0012  and  12%  thick  Joukowski 
airfoils  has  been  evaluated  with  global  PNS 
relaxation.  Analytic  or  conformal  mapping  [13] 
is  used  to  generate  the  requisite  metric  functions 
for  the  system  (1),  Figure  8.  Solutions  have 
been  obtained  for  fully  laminar  conditions  for 
Re  »  2000  to  7500.  Recirculation  is  evident  for 
the  Joukowski  airfoils  for  Re  >  2000.  For  the 
NACA  0012  configuration,  separation  is  not  evident 
for  Re  <_  5000.  Typical  laminar  solutions  are 
shown  in  Figures  (9a,  9b) .  The  laminar  stagnation 
point  results  are  also  in  close  agreement  with  the 
familiar  Navier-Stokes  (boundary  layer)  values  [4] . 

Triple  deck  analyses  have  recently  been  pre¬ 
sented  for  separation  on  cusped  and  sharp  trailing 
edge  airfoils  [14,  15] .  Estimates  of  incipient 
separation  as  a  function  of  Re  are  in  qualitative 
agreement  with  the  present  numerical  solutions, 
and  the  flow  behavior  near  the  wedge-like  trailing 
edge  is  also  reasonable,  see  Figure  10.  Further 
comparisons  are  given  in  [4] . 

Finally,  for  Re  »  5  x  105,  transition  to 
turbulent  flow  conditions  is  assumed  at  x/c  -0.32. 
The  two  layer  eddy  viscosity  model  should  be 
representative  of  the  turbulent  flow  behavior  as 
separation  does  not  occur  even  for  this  very  large 
value  of  Re  [16].  Comparisons  with  experimental 
results  and  earlier  calculations  [16]  are  quite 
reasonable,  see  reference  [4]  . 


Fig.  5b.  Trough  Pressure  Solutions  PNS  Solver: 
e » -0.03. 
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A  global  or  relaxation  procedure  for  the  PNS 
system  of  equations  has  been  developed.  First  and 
second-order  accurate  formulations  have  been 
presented.  In  the  latter  case,  a  staggered  grid 
is  considered  and  the  unknown  pressure  is  evaluated 
one  grid  point  upstream  of  the  velocity.  In  the 
former  case,  a  forward  pressure  difference  is 
implied.  The  effective  difference  forms  of  the 
vorticity  transport  and  poisson  pressure  equations 
have  been  derived  and  results  of  global  stability 
and  convergence  analyses  have  been  reported. 

Solutions  have  been  obtained  for  laminar  and 
turbulent  flows  where  strong  pressure  interaction 
and/or  axial  flow  separation  occurs.  The  full 
elliptic  pressure  interaction  is  accurately 
evaluated  and  with  the  local  pressure  interaction 
there  is  no  separation  singularity.  Procedures 
for  increasing  convergence  rates  have  been 
examined,  e.g. ,  multi-grid;  however,  further 
analysis  is  still  necessary,  see  [4]. 

Significantly,  the  differential  form  of  the 
poisson  pressure  equation  is  not  required  expli¬ 
citly  and  the  local  continuity  equation  is 
satisfied  exactly  at  all  points. 
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